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ABSTRACT 

Context. Planets are thought to eventually form from the mostly gaseous (~ 99% of the mass) disks around young stars. The density 
structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, far-ultraviolet (FUV), 
and X-ray wavelengths, as well as by the dust properties. 

Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the 
basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C + , O, H 2 0, and Ne + observed with facilities such as 
Spitzer and Herschel. 

Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and 
surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated 
using Pro toplanetary Disk Mo del (ProDiMo), with the recent addition of X-rays. 

Results. The abundance structure of neutral oxygen is very stable to changes in the X-ray and FUV luminosity, and the emission lines 
will thus be useful tracers of the disk mass and temperature. The C + abundance distribution is sensitive to both X-rays and FUV. The 
radial column density profile shows two peaks, one at the inner rim and a second one at a radius r = 5 — 10 AU. Ne + and other heavy 
elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from 
the lowest (Lx = 10 29 erg s~') to the highest considered X-ray flux (Lx = 10 32 erg s~'). FUV confines the Ne + ionized region to areas 
closer to the star at low X-ray luminosities (L x = 10 29 erg s~'). H 2 abundances are enhanced by X-rays due to higher temperatures in 
the inner disk than in the FUV only case, thus leading to a more efficient neutral-neutral formation channel. Also, the higher ionization 
fraction provides an ion-molecule route in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus 
providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. 

Key words, protoplanetary disks: X-rays - protoplanetary disks: disk structure, chemical composition 



1. Introduction 

New observing facilities in the past decade pushed our under- 
standing of protoplanetary disks from a rough picture of a ver- 
tically layered structure to a wealth of details on the composi- 
tion and two-dimensional structure of the gas and dust of these 
disks. In the infrared, the Spitzer Space Telescope performed 
systematic studies of nearby star-forming regions. The spectral 
energy distributions (SEDs) revealed the physical structure of 
disks, e.g., the presence of gaps, source-to-source variations and 
important statistics on SED types, which allow evolutionary sce- 

20 1 0) l . Extensive gas-phase 



narios to be built (cf. IMerin et al. 



emission line studies with Spitzer provide a first i ndication of 
chemical diversity across stellar spectral types (cf. Najita et al. 
2009] |Pontoppidan et al.|2010| |Lahuis et al.|2007[ ). In the near- 
infared, there have been large ground-based studies with, e.g., 
the VLT and the Keck telescope resolving the lines to study the 
kinematics (spatial origin of the lines in hot disk surfaces and 
winds) and their excitation mechanisms (thermal and fluores- 
cence; [Brittain et al. 2010, Pontoppidan et al. 201 1 Fedele et al. 
201 1} . In the past two years, the Herschel Space Observatory ex- 
tended the spectral window to the far-IR and submm, with spec- 



tral line scans within the Dust, Ice, and Gas in Time (DIGIT) key 
program (PI: N. Evans) for many Herbig disks, e.g., HD 100549 
( Sturm et al.|2010|l, water studie s of a few selected targets, such 
as TW Hya ([Hogerheijde et al. 201 1 J > from the Water in Star- 
forming Regions with Herschel (WISH) key program (PI: E.van 
Dishoeck), and large statistical gas surveys targeting the domi 
nant cooling lines (e.g.,|Mathews et aL2010; Meeus et al.|2012 



12; 

L2) 



|Dent & the GASPS team|2012||Riviere-Marichalar et al.|2012] 
in the Gas in Protoplanetary Systems (GASPS) key program (PI: 
B. Dent). The HIFI and some PACS line detections deal with the 
cold water, while the 63.3 //m H2O line possibly probes the inner 
water reservoir. 

Observational studies report trends in emission line strength 
with the irradiation of the central star. Gti del etaL] ( |2010| l find 
that the [Nen] emission at 12.81 yum correlates with the X-ray 
luminosity of the central star; the slope of the correlation for 

pTTo 



0.18. Pinte et al. 



and Meeus 



nonjet sources is 0.44 j_ 

|et al.| ( |2012| l show that the [0 1] 63 /zm line flux increases with 
stellar luminosity and that this is largely driven by the strength 
of the far-ultraviolet (FUV) luminosity. Riviere-Marichalar et al. 
( 20 1 2 ) report a tentative trend of the 63.3 fim water line flux with 
X-ray luminosity. 
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It is obvious from previous works that the most relevant pa- 
rameter for the SED is the total stellar luminosity; in most cases, 
Lx and Lfuv present only a small fraction of this and hence by 
themselves do not cause significant SED changes. However, re- 
cent near-IR high-contrast imaging with HiCIAO on the Subaru 
8.2 m telescope allowed probing of the inner disk structure on 
scales below 0.1" ( jThalmann et al.|2010||Hashimoto et al.|201 1) . 
This offers the possibility of a direct probe, e.g., the height of the 
inner rim as predicted by protoplanetary disk models. Thi et al. 
(201 1 ) predicted that the height of the inner rim could be a fac- 
tor two higher when the vertical hydrostatic structure takes into 
account the gas temperature at the inner rim. |Aresu et al. ( 201 l| l 
show that the height of the inner rim increases in disk models 
with the impinging stellar Lx- 

In the last decade, much theoretical progress has taken place 
to support the interpretation of the wealth of new observa- 
tions, specifically the gas observations. The foundation was laid 
by |Chiang & Goldreich] ( fT997] l, |D'Alessio efaX] (Q998), and 
Dull emond et al.| ( |2001| l, whose layered dust disk models have 
helped to interpret SEDs for irradiated disks (see Dullemond et 
al. |2007| for a review of disk structure models). Based on this, 
several groups focussed on the effects of FUV and X-rays on the 
thermal and chemical structure of the gas in a pre-prescribed pro- 
toplanetary disk model, e.g., [ van Zadelhoff et al.|(|2003|l, kamp 
& Dullemondl ( [2004l ), |Jonkheid et al.|(|2004|>, andlGlassgold et"a ~ 



0.01 



(|2004|l. In a subseq uent step,|Nomura & Millar| ( |2005| l, |Gorti & 
|Hollenbach| ( |2"00"8"l >, and |Woitke et al.| ( |2009| rstudied the chemi- 
cal structure of disks, while self-consistently solving for the ver- 
tical hydrostatic equilibrium using the gas temperature. Most re- 
cently, Aresu et al. ( 201 l| l performed an exploratory study to 
assess the relative importance of FUV and X-rays by expand- 
ing the existing Pro toplanetary Disk Mo del (ProDiMo) code to 
include X-ray processes. The modeling efforts that solve for the 
vertical disk structure are computationally expensive because the 
chemical networks, heating/cooling balance, 2D continuum ra- 
diative transfer, and hydrostatic equilibrium have to be solved it- 
eratively. Hence, these studies have largely focussed on a single 
representative disk model, or at most a handful of models vary- 
ing one specific parameter. However, large-model grids are in 
principle required to understand the potential diagnostic power 
of certain observables. So far examples for this are large grids 
(of the order of 200000-300000 model s) of parametrized dust 
disk models to study SED diagnostics (Robitaille et al. 2006) 
and of parametrized gas/dust disks to study the gas emission line 
diagnostics ( |Woitke et al.|201 l||Kamp et al.|201 1| >. 

In this paper, we perform for the first time an extensive anal- 
ysis of a grid of 240 self-consistent disk models (including the 
vertical disk structure) to study the effects of X-rays, FUV, and 
the relative importance of grain size and gas surface density 
distribution on the thermal, chemical, and physical structure of 
disks around T Tauri stars. This paper builds on the implemen- 
tation of X-rays into the ProDiMo code as described in | Aresu 



et al. (2011 1. While we discuss here foremost the physical and 
chemical structure of the models and how they change with ir- 
radiation, a companion paper ( Aresu et aL||2012| paper II), dis- 
cusses the power of line diagnostic in disentangling these effects. 
The paper is structured in the following way: Sect. [2] describes 
the updates on the code and the range of parameters used in the 
grid. The effects on disk temperature and density structure will 
be discussed in Sect. [3] The resulting distribution of various key 
species and its key reactions are extensively described in Sect. [4] 
Sect. [3] shows the radial column density distributions. Those are 
key in understanding the line profiles, which is the topic of the 
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Fig. 1. Input spectrum for model with Lx = 10 30 erg s" 1 and 
L FUV = 10 31 ergs-'. 



accompanying paper EL Sect.|6]summarizes the conclusions and 
implications of the paper. 



2. Updates on ProDiMo and the calculated grid 

The original ProDiMo code (|Woitke et aT]|2009|) i s based on 
the code as discussed by Kamp & Dullemond (2004). The orig- 



inal model includes (1) frequency-dependent two-dimensional 
dust continuum radiative transfer, (2) kinetic gas-phase and FUV 
photo-chemistry, (3) ice formation, and (4) detailed non-local 
thermodynamic equilibrium (non-LTE) heating and cooling with 
(5) a consistent calculation of the hydrostatic disk structure. 
The models are characterized by a high degree of consistency 
among the various physical, chemical, and radiative processes, 
since the mutual feedbacks are solved iteratively. |Aresu et al.| 
(2011 1 included X-ray heating and chemical processes, guided 
by work from |Maloney et al.| j!996| l, |Glassgold et al.| ( 2004| i, 
and |Meijerink & Sp aans (2005), while most recently the X-ray 



chemistry processes are updated following |Adamkovics et al. 
( |201 1) . This includes an extension of our chemical network with 
species such as Ne, Ar, and their singly and doubly ionized 
states, as well as other heavy elements. We added Ne + , Ne 2+ , 
Ar + and Ar 2+ in the heating/cooling balance and also imple- 



mented an extended sulfur chemistry (following Leen & Graff 

SO, SQ 2 , HS, and their re- 



1988, Meijerink et al. 



2008 



e.g.: 



lated reactions were added) to achieve a proper calculation of 
the sulfur-based species abundances. This allowed us to make 
more realistic predictions of, e.g., the sulfur fine-structure lines. 
A detailed description on the treatment of X-ray chemistry in the 
code is given in Appendix [A] 

Our input spectrum is composed of the stellar spectrum and 



a FUV excess as described in Woitke et al. (2009]), along with 
an X-ray component as outlined in |Aresu et al.| ( |2011) . The 
FUV excess (from 92.5 to 205 nm) spectrum is a power law 
with I V (A) cc A 12 . The X-ray spectrum (0.1 - 10 keV) has been 
chosen following Glassgold et al.| ( |2007| ) and is the same we 



used in Aresu et a 



( 201 \\ . It is a bremsstrahlung spectrum with 



I V (E) oc 1/E x exp(-£ IkT) (see Fig.[T|. We varied those quan 
tities that could potentially affect the penetration of X-ray and 
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FUV radiation and change the energy deposition through the 
disk. For example, we use two different values for the dust size 
distribution power law. The first value, generally used in the lit- 
erature for the ISM, is a power law index of a pow = 3.5. It is 
well known, though, that dust coagulation is the seed process 
that in timescales of a few million years leads to the formation 
of planetesimals. To model a disk that could have undergone 
dust coagulation, we also adopt a shallower power law index 
of fl pow = 2.5. For the same reason, we vary the minimum dust 
size: 0.1, 0.3, and 1 fim. These two parameters (a m [ n and a pow ) 
determine the FUV opacity and, as a result, impact the amount 
of energy reprocessed by the photoelectric heating. Ultimately, 
the parameters affect the gas thermal balance in the disk. The 
changing opacity and optical depth also impact the dust temper- 

atu re r dust . 

Aresu et al. ( 201 1) performed an exploratory study on the in- 



Table 1. Parameters used in the models. 



fluence of X-rays and FUV irradiation on the structure of proto 
planetary disks, for X-ray luminosities ranging from L x = 10 29 - 
10 32 erg s _1 . The range of values for this particular parameter is 
based on the outcome of the Taurus survey ( jGiidel et al.|2007 1, 
which showed that X-ray-active T Tauri stars emit in the range 
L x = 10 29 - 10 31 ergs" 1 . The model with L x = 10 32 ergs" 1 , 
on the other hand, is an extreme case with interesting implica- 
tions in terms of modeling. FUV excess emission from classic T 
Tauri stars between 6 and 13.6 eV is believed to arise from accre- 
tion spots on the surface of the star and stellar activity. The disk 
is truncated to several stellar radii from the star magnetic field, 
which channels matter toward the stellar photosphere. Shocks in 
the outer layers of the star cause the temperature to rise and emit 
FUV photons, which results in a FUV excess with respect to the 
stellar luminosity in the same band (Gullbri ng et al.|2000 1. Gorti 
& Hollenbach (2008 1 infer the FUV excess luminosity (91.2 
nm < hv < 205 nm) to be between Lpuv/L* ~ 10~ 2 - 10~ 3 . 
In our grid, we decided to scale the FUV flux so that we ob- 
tain the same range of luminosities as we use for the X-rays. 
The adopted values of Lpuv/L* are then 2.6xl0~ 2 , 2.6xl0~\ 
2.6xl0~ 4 , and 2.6xl0~ 5 . This way we can directly compare the 
input energy radiation in the FUV and X-ray band as well as 
its effects on the disk. We also explore two different values for 
the surface density power law distribution e, which is defined as 
£(r) = S(r )x(r/r )" £ : |Hayashi1 ( |1981| > derived the value e = 1.5 
in his model for the minimum mass solar nebular (MMSN) and 



their diagnostic value for our own solar system, while Hartmann 
( 1998 ) suggested e — 1 for objects older than 1 Myr. 

In summary, the current paper discusses the effects of vari- 
ations in the following parameters: X-ray luminosity (five val- 
ues), FUV luminosity (four), surface density profile (two), and 
dust size distribution, varying both the minimum grain size a m ,-„ 
(three) as well as power law indices (two). This yields a total of 
240 models; the summary of the model parameters is given in 
Table [1] 

3. Disk thermal and density structure 

The coupling of gas and dust to the radiation field determines the 
heating and cooling rates, which in return determine the pressure 
balance and structure of a disk. FUV and X-rays couple differ- 
ently to the gas: FUV is absorbed by dust grains and ejects an 
electron with a small surplus of energy in the form of kinetic en- 
ergy, which then consequently heats the gas. X-rays, on the other 
hand, are absorbed by the K-, L-, or M-shell of atoms, where no 
distinction is made as to whether this atom is part of a molecule, 
dust grain, or PAH. This assumption might overestimate the X- 
ray absorption cross section by a factor of approximately ~ 2 at 
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Fig. 2. Average spectral energy distributions of models with 
L x = (red), 10 29 (blue), 10 30 (orange), 10 31 (green), and 
10 32 erg s (yellow). 



energies E < 1 keV, since we do not include self-shielding ef- 
fects by large grains (see, e.g., Fig. 1 of Firema n|1974| l. On the 
other hand, not much is known about gas and dust phase elemen- 
tal abundances in disks. The observed abundance of, e.g., neon 
(dominating the X-ray absorption cross section between ~ 0.9 to 
1 .2 keV, and, as such, the total absorption of our 1 keV thermal 
source), already varies over a larger range (see Glassgold et al. 
2004, for a discussion on this topic). The electron that is ejected 
has a large energy capable of exciting, ionizing and heating the 
gas. The efficiency for X-ray heating is much larger, about 10 
to 40 percent, depending on the ionization fraction of the gas, 
while that for FUV heating does not exceed one to three percent. 
The coupling of X-rays to gas is weaker than for FUV, due to 
the smaller cross section for the absorption, but the larger effi- 
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Fig. 3. Total density: FUV luminosity increasing from Lfuv = 10 29 (left) to 10 32 erg/s (right). X-ray luminosity increasing from 
L x = (top) to 10 32 erg/s (bottom). 



ciency makes sure that it is an important heating source, even 
when the FUV excess is large. Their relative contribution to the 
energy budget will result in a different structure and emitted line 
fluxes. X-rays do not much affect the continuum fluxes, certainly 
less than a m j n , a pow , and e (see Fig. |2ji. It shows that the fluxes 
change at most by a factor of three between A = 5 - 100 yum due 
to X-rays. The models with Lx = 10 32 erg s are affected most, 
and these high X-ray fluxes are rarely observed for T Tauri stars. 



Woitke et al.|(|2009[) have already pointed out that the cooling 



cal structure of the disk. They present in their Fig. 9 the density 
structures for (1) a model where the gas temperature is decou- 
pled from the dust and (2) where the gas and dust temperature 
are coupled. Model 1 shows a very complex structure, where the 
gas is puffed up at the inner rim and a second bump in the den- 
sity is seen at a radius r ~ 10 AU. Aresu et aL"1 (201 1) extended 
this model by including X-rays and showed for the particular 
model of Woit ke et al.| (|2009) that the vertical extension inward 
of 10 AU becomes smoothed out for increasing X-ray fluxes and 



and heating balance is very important in determining the physi- 
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Fig. 4. Gas temperature distribution: FUV luminosity increasing from Lfuv = 10 29 (left) to 10 32 erg/s (right). X-ray luminosity 
increasing from L x = (top) to 10 32 erg/s (bottom). White contours represent the A v = 1 and 10. The black contours represent the 
gas temperatures. 

merges with the inner rim for the largest X-ray fluxes. A mini- always show the models with a m ; n = 0.1 jum, e = 1.5, and a pow 
mum in the vertical extension remains visible, though. 3.5. 



In Figs.[3]and|4] the density and temperature of the disk are 3.1. Density distribution 
shown for FUV luminosities ranging from Lfuv = 10 29 , 10 3 , 

10 31 , and 10 32 erg/s (left to right) and for X-ray luminosities Because density and temperature structure are directly related 

L x = 0, 10 29 , 10 30 , 10 31 , and 10 32 erg/s. There are 12 different and we assume hydrostatic equilibrium, many of the effects that 

possible variations with these X-ray and FUV luminosities (due are seen in the density structure will show up in the temperature 

to the variations in a m j n , a max , and e); as a baseline model, we structure as well. The FUV-only models show that the inner rim 
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Fig. 5. Scale height of the disk at z/r = 0.5. The black lines show the 12 different models (a m j n , a max , and e), at a particular Lfuv 
and Lx- The vertical green lines indicate a maximum or minimum in the scale height. The red line is the average of the 12 models. 



becomes higher and broader for increasing luminosity. The disk 
right behind the inner wall is shielded, and a second, vertically 
more extended density structure thus appears at radii between 
r — 2 and 6 to 30 AU (depending on the FUV luminosity). 

The X-rays have a different effect on the density structure. 
While the peak of the second bump in the density structure ap- 
pears at larger radii for higher FUV luminosities, the X-rays 
mostly affect the region within 5 AU. The second bump behind 
the inner rim expands to smaller and smaller radii for increas- 
ing X-ray luminosities (while keeping FUV the same), and the 
absolute densities become higher (from n ~ 10 6 to 10 s cm -3 ). 
Eventually, the puffed-up inner rim and the second bump merge. 
At the highest FUV luminosity, the merging of the two density 
structures occurs at the smallest X-ray luminosity. 

3.2. Temperature structure 

When only FUV is included in the models, we find that the tem- 
perature structure of the disk shows an inversion in the vertical 
direction: The temperature reaches a maximum at a certain rela- 
tive height z/r and shows a fast drop above it. For example, the 
models with L FUV = 10 29 erg/s show a region with temperatures 
T > 1000 K out to a radius r ~ 3 AU and between relative 
heights z/r ~ 0. 1 to 0.4. Noticeable is also that the hot region 
extends to a slightly higher relative height at the inner rim (up to 
z/r - 0.45) and at a radius r — 3 AU (out to z/r = 0.55). This is 
the result of shielding; the disk surface is less heated right behind 



the inner rim. For larger FUV luminosities, shielding effects are 
even more prominent and the drop in temperature right behind 
the inner rim is increasingly pronounced. 

Adding X-ray heating changes this picture. These models 
show a much more extended region with temperatures higher 
than T > 1000 K, even when the X-ray luminosity is much 
smaller than the FUV luminosity. At radii r < 5 - 10 AU, the 
inversion layer disappears and the temperature just smoothly in- 
creases toward larger z/r. At larger radii (r > 10 - 20 AU), 
the inversion layer only disappears when the X-ray luminos- 
ity is at least ten percent of the FUV luminosity. In models 
where the X-ray luminosity is much smaller, temperature dis- 
tribution is also affected (or even dominated) by the FUV. An 
example is Lfuv = 10 32 erg/s in combination with Lx = 10 29 
erg/s. Although the inversion layer disappears at small radii 
(r < 5 AU), it is still present at larger radii. 

3.3. Heating and cooling processes 

Temperature and density structure are in the end determined by 
the balance between the heating and cooling processes. Line 
cooling processes (except Lyman a and [OI] 6300A line cooling) 
are treated in an escape probability approximation. Before solv- 
ing the equations for statistical equilibrium, a global continuum 
radiation transfer calculation is done to estimate the background 
mean intensities for the radiative excitation and de-excitation 
rates. Other heating and cooling processes are approximated by 
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an analytical expression. Over 50 different heating and 50 dif- 
ferent cooling processes are included in the code. Except for 
the treatment of X-ray related heating processes, which are de- 
scribed in |Aresu efaX] ( |201 l| i, we refer to |Woifke et al.| ((2009 
201 1\ for a more detailed description. The reason that so many 
processes are included is that they each play a significant role, 
depending on the ambient density and radiation field (FUV, X- 
rays, cosmic rays), which vary over many orders of magnitude. It 
is beyond the scope of this paper to describe all these processes 
in detail for the models. However, we do present the dominant 
heating and cooling processes in Figs. B.l and B.2 for the mod- 



els for which we are showing the temperature and density struc- 
ture. 

First, we consider the FUV-only models. At the lowest FUV 
luminosity Lfuv = 10 29 erg s , there are three main heating 
process in the unattenuated part of the disk, namely, background 
heating by C + at the lowest densities, n H = 10 5 cirT 3 (dark 
blue), PAH heating within the inner rim and the second exten- 
sion of the disk (orange), and heating by carbon ionization in 
the outer disk (blue-green). Increasing the FUV field makes the 
picture more complicated. Heating by collisional de-excitation 
of H2 becomes important in the second bump (light blue), while 
background heating by Fell is dominant at the inner rim. Locally, 
other processes, such as photo-electric heating, are important. 
Going to the more shielded regions of the disk, there is roughly 
a three-layered structure of dominating heating processes: (1) in- 
frared background heating CO ro-vibrational transitions (black), 
(2) slightly deeper in the disk heating by thermal accomodation 
on grains (white), and (3) in the mid-plane cosmic ray heating 
(red). This three-layered structure becomes more confined to- 
ward the mid-plane for higher FUV luminosities. 

X-ray Coulomb heating dominates in the upper part of the 
disk, when the FUV to X-ray luminosity Lp^y/Lx < 1. When 
it is higher than one, X-ray Coulomb heating is restricted to the 
higher regions of the disk. Adding X-rays increases the temper- 
ature significantly in the disk, making [Fell] also dominant in 
large parts of the second-density extension at the lowest X-ray 
luminosity, Lx = 10 29 erg s _1 . For large X-ray luminosities, the 
structure of heating sources becomes less complicated. The en- 
tire unshielded part of the disk is dominated by X-ray Coulomb 
heating, followed by the three-layered structure described above 
(CO ro-vib, thermal accommodiation on grains, cosmic ray). For 
the highest two X-ray fluxes, an additional layer with X-ray H2 
dissociation heating as the main heating source is located on top 
of this three-layered structure. 

In the FUV-only case, there are three main coolants in the 
unattenuated part of the disk: C + line cooling (yellow), [Fell] 
line cooling (blue-green), and [OI] line cooling. The size of the 
region, where [Fell] line cooling dominates expands for higher 
FUV luminosities, while reducing the size of region where [OI] 
line cooling dominates. When X-rays are added, the region 
where C + line cooling dominates is pushed to the outer part of 
the disk. Temperatures are much higher in the upper part of the 
disk, and as a result Lyman a cooling (black) dominates in in- 
creasingly larger regions of disk, when X-ray luminosities be- 
come larger. At the highest X-ray luminosities, we find a three- 
layered structure of Lyman a cooling, [Fell] line cooling, and 
[OI] line cooling. The shielded region of the disk shows a lay- 
ered structure. CO ro-vibrational (red) and H2O rotational cool- 
ing (green-blue) are on top. Closer to mid-plane, several smaller 
regions have their own major coolant, such as HCN (purple), and 
HNC (orange). 



3.4. Vertical scale height 

Another way to show the combined effects of X-rays and FUV 
on the disk structure is by comparing the scale vertical height. 
Our definition is the same as the one used by Woitke et al. 
(2009), and is approximately (i.e., assuming z •« r) given by 



(If- 



2rcr 



(1) 



where H is defined as p(z) ~ p(0) exp(z 2 /// 2 ) and cj is the 
isothermal sound speed. As mostly the disk atmosphere is af- 
fected by the FUV and X-ray irradiation, we plot H/r at the 
relative height zjr = 0.5 in Fig. [5] In each panel, we over- 
plot the scale height for all the different parameters in black 
(twelve models in total), and the red line is the average. Four 
blue lines are overplotted, indicating flaring index H/r oc r p with 
p = 0.15 - 0.45, where the lines are normalized to the vertical 
scale height at r = 100 AU. 

A number of things stand out in these plots: (i) The vertical 
scale height of the inner rim is mostly affected by the X-rays, 
and almost no response is seen for different FUV luminosities. 
The scale height is H/r = 0.15 without X-rays at the inner rim, 
and it increases to ~ 0.4 at L x = 10 32 erg/s. (ii) The second 
bump also shows up when plotting H/r and slowly moves out- 
ward for increasing X-ray luminosities, (iii) The flaring index p 
is 0.25 in the outer disk (r > 100 AU), similar to the 2/7 ex- 



ponent found by Chiang & Goldreich (T997JI. The vertical scale 
height in the inner disk is elevated, compared to the outer disk, 
(iv) The H Q /r « 0.1 and 0.18 at r Q = 100 AU and zjr = 0.1 and 
0.5, respectively. For a direct comparison to Chiang & Goldreich 
( 1997} , we measure Ho/r at ro = 10 AU, because the approxi- 
mation of the disk dust temperature being vertically isothermal 
does not hold in our models at the larger distances. We find a 
value Ho/r » 0.04, which is very similar to Chiang & Goldreich 
( |1997| l, after correcting for the stellar mass and their definition 
of the scale height, (v) The merging of the inner rim with the 
second bump is clearly visible, as changes in the vertical scale 
height in the inner disk (r = 1 — 10 AU) become more gradual 
for larger X-ray luminosities. 



4. The chemical balance in the disk 

The combined effects of X-rays and FUV irradiation on the disk 
are discussed in this section. The density structure is in pressure 
equilibrium as discussed in the previous section and is altered 
due to irradiation effects. The chemical rates are density depen- 
dent, and as such the pressure balance also influences the abun- 
dances of the species. 

4.1. Electron abundances 

The main electron donor in a FUV-only chemistry is atomic car- 
bon, because the incident FUV radiation field is not capable of 
ionizing atomic hydrogen, leaving cosmic rays as the only ion- 
ization source of atomic hydrogen. As a result, the maximum 
electron abundance will not be higher than x e ~ 10~ 4 - 10~ 3 , 
as can be seen in Figs. [6] and Fig. |B.3| The maximum elec- 
tron abundances occur where the disk is unshielded to the ra- 
diation source. The inner rim has significant ionization fractions 
(x e ~ 10~ 4 - 10~ 3 ) all the way down to the midplane of the disk, 
while further out the electron abundances are only this high at 
densities n < 10 7 cm 4 . The drop of the electron abundances be- 
low x e < 10~ 6 - 10~ 5 for higher densities nicely coincides with 
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Fig. 6. Electron abundance: FUV luminosity is fixed at 10 31 erg/s. X-ray luminosity increases from L x = (left) to 10 32 erg/s (right). 
The black dashed line indicates Ay = 1. The blue lines mark the ionization fractions x e = 10 2 and 10 _1 . 




Fig. 7. Abundance structure of atomic (top) and molecular (bottom) hydrogen. FUV and X-ray fluxes are the same as in Fig. § The 
blue contours indicate the gas temperature r eas = 300 and 1000 K. The black line indicates 7d US t = 100 K. 



the A v = 1. Closer to the midplane, when the radiation becomes 
more shielded and densities are higher, the electron abundance 
exhibits a fast drop. This is the result of a combination of a de- 
creasing ionization rate, ( oc «, and an increasing recombina- 
tion rate, k lec oc « 2 . The abundance drop becomes more gradual 
at larger radii, which is a result of lower ambient densities and 
thus lower recombination rates. The region where the transition 
occurs is moving slightly toward smaller relative heights, when 
FUV fluxes are higher. The change in the transition location is 
not much due to the sharp increase of the density toward the 
midplane (k rec oc « 2 ): A relatively small increase in the density 
already compensates for larger FUV fluxes. 

This picture changes quite significantly when the radiation 
fields also contain X-rays. The absorption of an X-ray photon by 
an arbitrary species creates a fast electron that can produce many 
ionizations, e.g., a 1 keV electron is able to produce ~ 27 hydro- 
gen ionizations. As a result, the electron abundances are much 
higher when X-rays are included. In the unattenuated part of the 
disk, ionization fractions are of the order x e ~ 10~ 2 - 10 (con- 
tours for x e = 10~ 2 and 10~' are indicated). They can even be 
higher than x e > 10~' at very small radii (r < 0.5 AU) and high 
relative heights {z/r > 0.5). Here, low densities («h < 10 7 cm -3 ) 
and high temperatures (T » 5000 K) reduce the recombination 
rates. Recombination rates decrease with temperature up to T ~ 
10000 K because they are dominated by radiative recombination, 
which scales as k KC oc n e n-,(T I \0 A y Xa,i , with X la( j ~ 0.6 - 0.9. 



We find slightly higher ionization fractions than, e.g., Glassgold 
et al.| ( |2QQ4] >. Their calculation at r = 1 AU (see their Fig. 4) is 



truncated at vertical col umn density Nu ~ 2 x 10 18 cm 2 and 



density n H ~ 10 cm . Ercolano etal. (20081 (their Figs. 1 and 



2) show calculations for lower column densities and densities, 
and they find ionization fractions slightly higher than x e > 1 at 
0.07 AU at vertical column densities Nu < 10 16 cm" 2 . Both cal- 
culations are consistent with those presented here. The electron 
abundance structure looks a little counterintuitive, especially for 
the lowest two X-ray fluxes L x - 10 29 and 10 30 erg/s. The region 
with very high (x e « 10~ 2 - 10~') becomes smaller for increas- 
ing FUV luminosities (Fig. |B.3[ >. This is because the inner rim is 
puffed up more and even merged with the second bump for the 
highest FUV fluxes. The densities are elevated to «h = 10 8 cm -3 
in these regions, and the recombination rates are orders of mag- 
nitude higher as a result. The outer disk is shielded from radia- 
tion. The regions with very high ionization fractions (x e > 10~ 2 ) 
are of similar size for all FUV luminosities only for the highest 
X-ray luminosity (see Fig. B.3 i. 



When the ionization fraction is increased by orders of mag- 
nitude, the gas shows a radically different chemistry, as it will be 
dominated by ion-molecule reactions. This will be discussed for 
some of the key species, and water in particular. 
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Fig. 8. Abundance structure of H + , H , HI,, and H3. FUV and X-ray fluxes are the same as in Fig. ^ Black contour indicates 
n H — 10 5 crrT 3 . 



4.2. H and H 2 abundances 

The abundance structure of atomic hydrogen is shown in the top 
panel of Fig. [7] (and also in Figs. |B.4| and |B.5| l. The highest 
abundances (x H > 10~ 2 ) are obtained in the unshielded region 
of the disk, where FUV and X-rays have high fluxes. At radii 
r < 0.6 AU, the disk is atomic all the way down to the mid- 
plane of the disk, as it is directly exposed to the central source. 
At larger radii, the atomic fraction drops below xu ~ 10~ 2 and 
makes a transition to molecular hydrogen at a relative height of 
z/r ~ 0.1 between r ~ 0.6 - 5 AU. The transition occurs at in- 
creasingly larger relative heights at larger distances to the central 
source, and it is at approximately z/r ~ 0.4 - 0.5 (depending on 
the X-ray flux) at a distance r ~ 100 AU. The transition occurs 
at slightly smaller relative heights at larger X-ray fluxes. 

The formation of H2 on grains is extremely efficient (~ 1) 
when the dust temperatures are not high (7d USt ~ 10 - 50 K 
( |Cazaux & Tielens] 2004 1. At dust temperatures higher than 
Tdust > 100 K, the H2 efficiency drops rapidly and H2 forma- 
tion on dust does not occur at Td ust > 1000 K. H2 is also formed 
through the H route, H + H — > H2 + e , especially when X- 



rays are present, and the fractional abundance of electrons is rel- 
atively high (and thus also H ). It is possible to maintain an effi- 
cient route to form molecular hydrogen in the gas phase at high 
temperatures, T gas > 300 K and rd ust > 100 K (see contours 
in Fig. IT}. Also FUV is much more efficient in destroying H2, 
than X-rays. The FUV photo-dissociates H2 (H2+ FUV photon 
— > Hj — > H + H). X-rays predominantly ionize molecular hydro- 
gen indirectly by collisions with fast electrons produced after an 
X-ray absorption (H2 + e~* — > + e~ + e~*), while only a small 
fraction of Fb is dissociated in this process. After ionization, H2 
is able to reform H2 through Hj +H — > H2 +H + , but it will also be 
able to form H3 through the reaction Hj +H2 — > Ht +H. This par- 
ticular species is key in forming molecules through ion-molecule 
reactions. Because X-rays are destroying H2 less efficiently and 
can provide an environment to form H2 in the gas-phase, a fast 
ion-neutral chemistry results, which makes it possible to form 
molecules (and maintain significant abundances) at much higher 
temperatures than when only FUV is present. 

The bottom panel of Fig.|7]presents the abundances of H2. It 
shows that the H2 abundance in the warm atmosphere of the disk 
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is largest for the highest X-ray fluxes (as high as jch, ~ 10 in 
the region directly exposed to the central source at r — 0.6 AU). 
This is the result not only of the processes described above, but 
also of the scale height of the inner rim, which increases with 
X-ray luminosities, an effect described in the previous section. 



4.3. HP, H~, H+, and H+ abundances 

The main drivers of the chemistry in the atmospheres of disks 
exposed to irradiation are ionic species, starting with H + , HJ, 
and H3 (see also the discussion of the water chemistry). The H 
abundance structure is shown for completeness, since it can be 
important in the production of H2 (see Fig. [8] and Figs. B.6 



to 



The abundance structure of ionized hydrogen shows abun- 
dances as high as xh+ ~ 10~ 3 - 10~ 2 , even if there is only FUV 
irradiation. FUV is not able to photo-ionize atomic hydrogen, 
and hydrogen ionization is solely done by cosmic rays. This only 
happens at very low densities (n < 10 4 cm 4 ), though. Once 
densities of the order of n ~ 10 5 cm 3 are reached, the ionized 
hydrogen fractions are closer to x e ~ 10~ 5 - 10~ 4 or smaller. 
The ionization balance is regulated by the charge exchange re- 
action H + + <-> H + + O. The backward reaction has an energy 
barrier of T — 227 K, so the ratio of + with respect to H + de- 
creases with temperature. The charge exchange rates are much 
faster than the photo-ionization reaction or the recombination 
rates of the species. The amount of H + is thus directly coupled 
to the production rate of + , which is also produced by cosmic 
ray ionization. Even when X-ray ionization determines the ion- 
ization balance of the gas, the charge exchange rates are much 
faster than the primary and secondary X-ray ionization rates of 
the species, and the ion abundances are directly coupled. X-rays 
produce much higher fractions of H + (as high as x^+ > 0.1 at 
the lowest densities). The abundance structure strongly resem- 
bles that of the electrons, although the H + abundance drops faster 
closer to the mid-plane, an d other species such as N a + take over 



Adamkovics et al. 
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as electron donor (see also 

Once the gas gets more shielded to the radiation, it has higher 
abundances of molecular species (see, e.g., Fig. |B.5| >. In those 
regions, H + can also be produced by other reaction paths, such as 
CO + + H — > CO + H + when only FUV is present and secondary 
X-ray ionization of H2 when X-rays are present as well. Once 
cosmic ray ionization is the dominant source of ionization, H + 
is produced by cosmic ray ionization of H2 or charge exchange 
with He + . 

The H ion is produced by cosmic ray ionization (H2+CR — > 
H + + H ) or by radiative recombination (H + e — > H +hv). This 
last reaction is fairly slow at low temperatures due to temperature 
barriers, but it becomes important in the warm part of the atmo- 
sphere of the disk. When only FUV is irradiating the disk, the H 
abundance is high (xh — 10 -11 - 10 10 ) at the inner rim and in 
the second bump extending out of the disk, and to a lesser extent 
in the transition zone from atomic to molecular gas. There are 
two reasons that there is not as much H at high altitudes in the 
disk. The first one is that the electron densities are two to three 
orders of magnitude lower. The second one is that the ambient 
temperatures are between 50 and 100 K, which is not very favor- 
able as the reaction rate contains an energy barrier. At the inner 
rim and in the second bump, both temperatures (r > 1000 K) 
and densities (n > 10 7 cm' 3 , and thus also electron densities) 
are much higher. When X-ray irradiation plays a role, the abun- 
dances structure changes a lot. The temperatures are in excess of 
T > 1000 K to large relative heights (see Fig. [4}, and electron 



abundances are much higher (due to higher ionization rates and 
slower recombination rates). Especially when X-ray luminosi- 
ties exceed Lx > 10 3() erg/s, the relative abundances of H~ are 
one to three orders of magnitude higher (close to x^- ~ 10~ 8 ). 
In the FUV case, the FF route to form H2 is at least three orders 
of magnitude lower than the formation route on dust grains. It 
becomes a significant contributor (~ 30%) to the total H2 forma- 
tion rate at the highest X-ray luminosities when these extreme 
(xu- ~ 10~ 8 ) abundances are reached, and at the same time, H2 
formation on dust is quenched by high dust temperatures. 

is formed by cosmic ray ionization (H2 + CR — > + e") 
in the regions where the gas is molecular and shielded from FUV 
and X-ray radiation. When the H2 abundance is high (xh 2 ~ 0.5), 
however, the HJ is quickly transformed to HJ (Hj + H2 — > H3 + 
H). It is above this layer, that the HJ abundance is increasing. 
Here, the gas is more exposed to radiation (FUV or X-rays). If 
X-rays are present, H2 is ionized through secondary ionizations. 
In the FUV-only case, where we do not have this last reaction, 
other routes contribute to the production of HJ . One example is 
S + +H2 — > Hj + S. In the unshielded regions of the disk, where H2 
abundances are low (xh, < 10~ 6 ), the only efficient way to form 
Hj is through radiative association H + +H — > H^ + hv. However, 
this is only efficient at relatively high temperatures (T > 5000 K) 
and significant abundance levels of H + . In Fig. [8] the cation H + 
is shown to be more abundant by orders of magnitude when X- 
rays are present. It is obvious that the production of is much 
more efficient as well, and hence abundances xm ~ 10~ 8 - 10~ 7 
occur in the high warm atmosphere of the disk. In the FUV-only 
case, there is only a significant abundance (xh+ > 10~ 12 ) in the 
inner rim and the second bump. 

The dominant way to form H3 is through the ion-molecule 
reaction H£ + H2 — > H^ + H. As a result, the H/H2 transition 
layer and the region below is very suitable to produce large abun- 
dances: H2 is present in reasonable amounts (x H , > 10~ 5 - 10" 3 , 
see Fig. [Tji, and radiation is available to produce Hj (see previ- 
ous paragraph, and Fig. [81. Once the HJ drops below an abun- 
dance xh+ < 10~ 12 , the H^ abundance also drops to values below 
xh+ < 10 -11 . Consequently, there is a distinct layer where HJ is 
formed in the most optimal way. We also saw that the H} is 
more efficiently formed when X-rays are abundantly present in 
the disk. The HJ abundances thus also significantly increase for 
larger X-ray luminosities and can be as high as xh+ ^ 10~ 8 . In 
the FUV-only case, the size of the region where H3 is present 
in large abundances is smaller, and the maximum abundance is 
at least an order of magnitude lower compared to the models 
that include X-rays. It is a key species in driving ion-molecule 
chemistry, and so the chemistry with and without X-rays will 
obviously yield significantly different abundance structures. 

4.4. C 2+ , C\ C and CO abundances 

The ionization potentials (IP) of the species C and C + are 
IP = 1 1.26 and 24.38 eV, respectively. While FUV photons can 
ionize neutral carbon, X-rays (or collisions with the resulting 
fast electrons) and/or cosmic rays are needed to reach a higher 
degree of ionization. As a result, the abundance patterns of the 
two species C + and C 2+ respond radically differently to the in- 
cident radiation field on the disk. The C 2+ abundance pattern 
is, of course, strongly correlated with the X-ray luminosity, and 
C 2+ is absent in models without X-rays. The region with sig- 
nificant C 2+ abundances (x c i+ > 10~ 9 - 10~ 8 ) extends to larger 
radii and smaller relative heights for higher X-ray luminosities, 
simply because there is more ionizing radiation available (and 
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Fig. 9. C 2+ , C + , C, CO abundances. The black dashed line indicates log(y/n H ) = -5. Fluxes are the same as Fig. [6] 



not inhibited by dust, since the dominant opacity is caused by 
the gas). However, the FUV plays an important role in shaping 
the disk and therefore indirectly affects the abundance distribu- 
tion as well. The bump in the total H number density (see Fig. 
[3]) extends to larger relative heights for larger FUV luminosi- 
ties, thereby increasing the density at larger relative heights. This 
causes a larger attenuating column between the central source 
and the outer part of the disk, thus reducing the amount of ioniz- 
ing photons. Furthermore, it allows H2 to form at smaller radii. 
H2 is very efficient in reducing C 2+ by charge exchange reac- 
tions, such as C 2+ + H 2 -» C + + H. As a result, the FU V lu- 
minosity confines the C 2+ to smaller radii (see Fig. B.lOi. C + is 
present throughout the unattenuated part of the disk. It has an 
abundance xc* ~ 10~ 4 down to the mid-plane in the inner rim 
(r < 0.6 AU). Beyond this radius, carbon is only significantly 
ionized above a relative height z/r ~ 0.1 (at r ~ 0.6 AU), in- 
creasing to z/r ~ 0.4 at r ~ 200 AU. The width Az/r, over which 
the transition from C + to C and CO occurs, is larger toward the 
outer regions of the disk. Absolute densities of C + and e~ in 
the inner disk are larger, and recombination rates scale with n 2 . 
As a result, the transition from ionized to neutral carbon occurs 



abruptly, similar to the abundance drop of electrons (at approx- 
imately z/r ~ 0.01 at r = 1 AU). However, in the outer disk 
(r = 200 AU), the transition stretches out from z/r ~ 0.04 to 
0.01. Although there are variations of a factor a few in the abso- 
lute abundance of C + , the overall appearance of the abundance 
structure is not significantly affected for different values of the 
X-ray and FUV luminosity (see Fig. |B.ll) . 

The neutral carbon abundance pattern is very much affected 
by variations in the X-ray and FUV luminosities. In the case 
where only FUV is irradiating the disk, one finds a clear transi- 
tion from C + to C to CO (see Fig.|9]>, which is expe cted in a FUV- 
dominated PDR (cf. |Hollenbach & Tielens|1999l >. To guide the 
eye, we added a contour with value log(x/nu) = -5 (with x the 
Draine field). This contour follows the C + /C/CO transition very 
well, indicating PDR physics. Neutral carbon is confined in a 
layer between the C + and CO. This picture is more complicated 
when X-rays are added. Perhaps not entirely intuitive, the over- 
all trend with increasing X-rays is that neutral carbon is more 
and more abundant in the unattenuated part of the atmosphere 
(as high as xq ~ 10~ 5 - 10~ 4 ). When X-rays are present, neu- 
tral carbon is higher by at least one to two orders of magnitude 
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Fig. 10. 2+ , + , and O abundances. Fluxes are the same as in Fig. ^ The black contour in the lower panel indicates the critical 
density of the [OI] 63 yum line. 



in abundance compared to models that exclude X-rays. The C 
and C + ionization rates are of course higher, but the recombina- 
tion rates are also higher owing to the increased electron abun- 
dances (although the higher temperatures reduce the recombina- 
tion rates, as already pointed out in Sect. 4.1 1. This is because 



electrons are now also produced by, e.g., ionization of hydrogen 
and helium. The fact that C + and C coexist in more or less equal 
amounts, when X-rays are dominating the ionization fraction of 



gas clouds, was already pointed out in papers by Maloney et al. 
1996 ) (see their Figs. 3a and 4a) and Meijerink & Spaans ( 2005 ) 



(see their Figs. 3 and 4). 

The bulk of the CO is situated in the shielded, optically thick 
part of the disk, and the total CO gas mass does not change 
when varying the FUV and X-ray radiation field. This molecule 
could thus serve as a tracer of the total mass of a protoplanetary 
disk, keeping in mind that the conversion factor depends on the 
amount of ice formation. A few details should be noted. There is 
a significant correlation between CO and OH in the unattenuated 
part of the disk (compare Fig. B.13 with B. 17 1. A route to form 



CO is through the reaction C+OH -> CO+H. While it is possible 
to obtain an abundance level of xco ~ 10~ 7 through this channel, 
it will never go up to xco ~ 10~ 4 , because the OH abundance is 
xoh ~ 10~ 7 . As OH becomes less abundant for higher FUV lu- 
minosities (when X-rays are fixed), the CO also becomes less 
abundant and slightly more confined to the midplane. In addi- 
tion, the the average abundance of CO drops slowly with increas- 
ing X-ray luminosity at the inner rim by destruction through ion- 
molecule reactions, such as CO+He + — > C + -i-O+He. This could 
have an effect on the ro-vibrational lines that are predominantly 
produced in these regions. 



4.5. 2+ , + , and O abundances 

The ionization potentials of O and + are IP = 13.68 and 
35.12 eV. Even the IP of neutral oxygen is above the threshold 
of IP - 13.6 eV, where neutral hydrogen blocks the radiation 
efficiently. The implications of this are seen in Fig. 10 where 



the abundances of neutral, singly and doubly ionized oxygen 
are shown for a FUV luminosity Lfuv = 10 31 ergs -1 and all 
considered X-ray luminosities except Lx - 10 32 ergs -1 . The 



second row of Fig. 10 (and also Fig. B.15 for variations with 
FUV luminosity) shows the abundance structure of singly ion- 
ized oxygen. Unlike the C + , the abundances of + do not ex- 
ceed xo+ ~ 10~ 8 - 10 in the FUV-only case. The abundance of 
ionized oxygen is mostly set by the charge exchange balance 
with atomic hydrogen, while the ionization of both hydrogen 
and oxygen is entirely due to cosmic ray ionization. When X- 
rays are added, the abundances increase three to four orders of 
magnitude compared to the FUV-only case, even at low X-ray 
luminosities (Lx = 10 29 erg s~'). Important to note is that the + 
is confined to smaller radii for larger FUV luminosities, whic h is 
also the case for C 2+ . This was discussed earlier in Sect. l4.4l and 
explained by the larger recombination rates due to higher densi- 
ties in the upper atmosphere of the disk and the lower ionization 
rates due to shielding. The confinement to smaller radii is even 
stronger for the case of 2+ , because this species is not affected 
by cosmic rays. 

The bottom panel of Fig. [TTJJshows the neutral oxygen abun- 
dance structure. The critical density is n cl ; t = 5 x 10 5 cm -3 
for the [OI] 63 /mi line, which is indicated by a black dashed 
contour. This species does not change with X-ray luminosity. 
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Fig. 11. OH, H2O, OH + , H20 + , and H30 + abundances. Fluxes are the same as in Fig. [6] The white contours indicate gas temperatures 
of 250 and 2000 K. The red contour indicates the number density «h = 10 6 cm 3 . 



Neutral oxygen is the dominant oxygen carrier throughout the 
disk, except for those regions, where large fractions of the gas 
are frozen onto dust grains (i.e., the midplane). This species is 
very insensitive to changes in FUV and the ambient chemistry 
(see Fig. |B.16| l. The bulk of atomic oxygen is in LTE for the 
commonly observed [OI] 63 and 145 /vm fine-structure lines and 
thus sensitive to the average temperature of the region that it is 
probing. As a result, the [OI] fine-structure lines have the poten- 
tial to probe the total energy budget of the gas (the combined 
FUV and X-ray luminosity irradiating the disk). There will be 



an elaborate analysis of the oxygen fine-structure line emission 
in paper II. 



4.6. OH, OH + , H 2 0, H 2 + , and H 3 + abundances 

The water chemistry and its related species are significantly af- 
fected by X-ray irradiation. The physical circumstances in the 
disk change with radius r and height z due to radiation shield- 
ing and large differences in chemistry, which means that the 
X-ray and FUV energy deposition per particle, Hx/n, and^/n 
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Fig. 12. Abundance structure of selected heavy elements: Ne + (top), Ne , Ar + , and Ar (bottom). The FUV luminosity is fixed at 
^fuv = 10 31 ergs -1 . X-ray luminosities vary between Lx = 10 29 and 10 32 ergs -1 . 



(with x — 1 the interstellar radiation field as defined by |Draine 
1978||Draine & Bertoldi|1996[ and integrated between 91.2 and 
205 nm), and their relative ratio change over orders of magnitude 
in the disk. As a result, the main chemical pathways forming the 
molecules OH and H2O change throughout the disk and the grid. 

X-rays heat the gas and therefore the reaction rates with ac- 
tivation barriers increase. They also drive the ionization, which 
enhances formation routes through ion-molecule reactions. The 
formation of the neutral species OH and H2O is either by the 
neutral-neutral reactions, H2 + O — > OH + H, and H2 + OH — > 
H2O+H, or by recombination of ionized species, such as H30 + + 
e" -> H 2 + H and H 2 + + e " -> OH + H. Which of these path- 
ways are dominating depends thus strongly on the temperature 
and ionization fraction. The routes through neutral-neutral reac- 
tions have activation barriers, making them efficient only when 
gas temperatures are sufficiently high, i.e., T > 200 - 300 K. 
One way to increase the neutral-neutral route without the need of 
higher temperatures is by lowering the activation barrier, which 
can be done by exciting H2 to a higher vibrational state after ab- 
sorption of a FUV photon. At high relative heights (z/ r > 0.5), 



this reaction will dominate when only FUV photons are present. 
The ion-molecule chain is started with cosmic ray or X-ray ion- 
ization of H2 and thus requires significant ionization rates. 

It turns out that the neutral-neutral reaction pathways dom- 
inate in the inner part of the disk (r < 5 AU). The X-rays heat 
regions deeper in the disk, resulting in the production of a thick 
warm water layer. Although there is a larger fraction of X-rays 
ionizing the gas than when only FUV photons are present, the 
X-ray heating causes the neutral-neutral reactions to dominate. 
It is remarkable to note that in the outer region (r > 20 AU) of 
the disk the ion-neutral reaction network dominates. This is be- 
cause temperatures are too low for the neutral-neutral reaction 
pathways to occur. These combined effects become apparent in 
the distribution of the water throughout the disk, as shown in the 
second panel from the top in Fig 



11 



(see also Fig. 



BT8] for a 
In 



more elaborate view on the effects of both FUV and X-rays) 
the models without any X-rays (left-hand side), the water shows 
a strong abundance peak at the inner rim (all the way down to the 
midplane of the disk) and a warm water layer at a relative height 
z/r m 0.1, with an abundance of xh,o ~ 10~ 7 , and a second wa- 
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ter reservoir higher up in the disk. This reservoir is located at a 
relative height z/r ~ 0.015 at small radii (r ~ 1 AU), and at a rel- 
ative height of z/r ~ 0.5 at a radius of r ~ 100 AU, thus follow- 
ing the flaring of the disk. These different water reservoirs were 
already noted by|Woitke et al. (2009 1 for Herbig AeBe stars. An 



elaborate discussion on the formation of warm water reservoirs 
through (the dominant) neutral-neutral reactions and formation 
on dust grains in the inner regions of an X-ray irradiated disk is 
also outlined by Glassgold et al. (2009). 

In the models where only FUV is included, the second layer 
is confined to the regions where temperatures are between T ~ 
300 - 1000 K. The regions at higher temperatures are exposed 
to FUV fluxes that essentially destroy the water faster than it 
is formed. When X-rays are included, the second water reser- 
voir extends to increasingly larger radii. Although there is water 
present in the FUV-only models out to radii r ~ 200 AU, the 
abundances are ultimately one to two orders of magnitude higher 
when X-rays are added. The additional water is increased since 
the ion-molecule chemistry is very effective in forming the wa- 
ter. This is illustrated in the third to fifth panel of Fig. 1 1 (and 
also Figs. B. 19 B.20 and B. 21 1, where the changing abundance 
structures of the ionic species OFF, H20 + , and H30 + are shown 
for various combinations of X-ray and FUV luminosities (ap- 
pendix A only). The ionic species extend to larger radii for larger 
X-ray fluxes, but larger FUV fluxes decrease the extent. The fact 
that FUV confines the ionic species to smaller radii is an indirect 
effect, because FUV tends to puff up the upper layers of the disk 
and the density becomes higher at larger relative heights (see 
Fig. [3]). These larger densities increase the recombination rates 
and reduce the abundances of these ionic species. 

The first layer described above is located in those regions of 
the disk, where the conventional routes do not work anymore. 
Temperatures are too low and the ionization fraction is small. 
This is the region where more exotic reactions take over, such as 
NH2 + NO — > N2 + H2O. This layer has a smaller vertical extent 
in relative height for larger FUV fluxes. On the other hand, the 
second separate layer becomes thicker when X-rays are added. 
When the Lx/Lfuv ^ 1 ratio, these two layers merge. The rea- 
son for the merging is the enhancem ent of the tra nsi ent spe cies 
OH\ H 2 + , and H 3 + (see also Fig.[Rl^[B50larid^2l). As 



mentioned earlier, these species react sensitively to the presence 
of an ionization source and become increasingly abundant with 
higher X-ray fluxes. 

The situation with OH is very similar. OH becomes more 
abundant in the outer disk, when X-rays are added. OH is located 
at higher altitudes in the disk (although there is also a large over- 
lap with the regions where the water is located). For this reason, 
the effect of the X-rays on the abundance in the outer disk is not 
as pronounced as for water. The H2 used in the neutral-neutral 
formation route of OH is on average at higher temperatures, so 
the relative contribution of the ion-molecule route is smaller. 

4.7. Ar + , Ar 2 ' 1 , Ne + , and Ne 2+ abundances 

These species have ionization potentials that are larger than 
that for atomic hydrogen, with IP = 21.56, 40.96, 15.76, and 
27.63 eV for Ne, Ne + , Ar, and Ar + , respectively. The only way 
to produce these species is through direct (photon absorption and 
Auger effect) or indirect (fast electron collisions) ionization by 
X-rays. The presence of these ionized species is thus a direct 
result of the X-ray irradiation of the disk. 

The models include a thermal source of temperature, T x = 
1 keV. As a result, the cross sections for neon are more favorable 
for direct ionizations than for argon. The cross section for the 



absorption is located at approximately ~ 1 keV for neon, while 
for argon it is ~ 4 keV. On the other hand, the rates for secondary 
ionizations are a little higher for argon, cr(A)/cr(H) =1.1, 0.48, 
3.7, and 1.8, with A equal to Ne, Ne + , Ar, Ar + , respectively. 
Consequently, the total ionization rates are very comparable. 
However, there is one big difference in the chemical network 
of these species. The charge transfer rate of Ne + with H2 is 
very low (k < 10~ 14 cm 3 s~'), while the other ions have sig- 
nificant rates for charge exchange with molecular hydrogen. The 
Ne + abundance structure thus extends to much smaller relative 
heights, down to z/r ~ 0.01, than Ar + , which is abundant above 
z/r ~ 0.01 - 0.06 at radii r > 10 - 200 AU. The abundance 
structure of Ar + is very similar to those for Ne 2+ and Ar 2+ . 

5. Radial column density profiles 

Although the abundance distributions of species well illustrate 
the change in the chemistry due to different combinations of X- 
ray and FUV luminosities, it is necessary to look at the integrated 
properties of the disk that are observed by our telescopes. The 
focus in this is on the species C + , O, H2O, and Ne + , as their line 
intensities and line profiles are extensively discussed the paper 
II. The figures show the average column density in yellow, while 
the black lines represent the 12 models at fixed FUV and X-ray 
luminosities. Changing parameters other than those for FUV and 
X-rays does not affect the results si gnifica ntly. 

C + column density profile (Fig. B.26): The column density 
of C + is strongly related to the inicident FUV flux at the inner 
rim, where it increases from Nq* ~ 10 17 to 10 19 crrT 2 , when 
the FUV luminosity increases from Lfuv = 10 29 to 10 32 ergs -1 . 
Right after the peak at the inner rim, there is a sharp drop in 
the column density, which is caused by the inner rim casting a 
shadow and therefore reducing ionizing radiation. A second peak 
in the column density profile is seen between r ~ 5 - 10 AU. 
The peaks moves to larger radii for higher FUV luminosities. 
The X-rays affect the column density profile in a different way. 
They reduce the minimum column density right after the inner 
rim and smooth the radial distribution. The combination of the 
highest FUV and X-ray fluxes gives the flattest distribution. 

O column density profile (Fig.\l3\: The column density pro- 
files of neutral oxygen do not change as much as those for C + . 
The maximum column density, No ~ 3 - lOx 10 19 cm~ 2 , is 
located right behind the inner rim and only varies with a factor 
of at most ten, which is again caused by the FUV irradiation. 
The width of the oxygen column density peak broadens a bit for 
higher FUV and X-ray fluxes. This makes the mininum in the 
column density less prominent in the radial column density pro- 
file. Because the column density is not very affected by radiation, 
it is a clean probe of the properties of a disk temperature ([OI] 
fine-structure lines), since there are no strong dependencies on 
uncertainties in the chemical network. 

H2O column density profile (Fig. \14\ : The highest water col- 
umn densities, N Hl o ~ 10 20 - 10 22 cm 2 , are found at the inner 
rim, r ~ 0.5 AU. As discussed, the water abundances are en- 
hanced by higher temperatures and higher ionization fractions 
throughout the disk. The FUV counteracts this to some extent, 
and the most favorable situation is a high X-ray to FUV luminos- 
ity ratio. FUV irradiation puffs up the inner rim, and shields the 
outer disk, thus allowing less ionizing radiation to penetrate into 
the outer regions of the disk. As a result, it confines the water to 
smaller regions of the disk. The bottom panel (Lx = 10 32 erg s _1 ) 
of Fig. 14 shows that the column density distribution is flat for 
the lowest FUV flux, dropping steadily as a function of radius 
for the highest FUV flux. The only region in the disk where the 
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Fig. 13. Radial column density distribution of O. FUV and X-ray fluxes are the same as Fig. [3] 



column density becomes smaller is at the inner rim due to re- 
actions with ions such as H2O + He + — > H + + He + OH. The 
water column density is almost two orders of magnitude smaller 
for the models with the highest X-ray luminosities compared to 
those with only FUV. 

Ne + column density profile (Fig. B.27): Ne + is only pro- 
duced by X-rays and therefore only the models that include X- 
rays show significant column densities. The profiles show a peak 
in the column density at the inner rim and a second bump at a 
few AU. The second bump smoothes out for larger X-ray fluxes, 
while the maximum column density increases from A^e + ~ 10 15 
to 10 17 cirT 2 . The FUV tends to confine the Ne + to smaller radii 
for the lower X-ray luminosities, which becomes apparent in the 
line profile (see paper II). 



6. Conclusions 

In this paper, we discussed the combined effects of FUV and X- 
rays on both the density structure and the thermal and chemical 
balance of disks around T Tauri stars for an expected range of 
parameters (dust size distribution, density profile, etc.), yielding 



a total of 240 models. Here we highlight the main results and a 
few implications: 



6.1. The disk thermal and chemical structure 

Temperature structure: The extent of the disk where the temper- 
ature is higher than T > 1000 K is much larger when X-rays 
are included. X-rays have a much higher heating efficiency than 
FUV, 30 - 50% compared to < 3%, respectively. 

Density structure: Increasing FUV luminosities does not 
change the scale height of the inner rim; it only alters the width 
and height of the second bump in the disk that is created at inter- 
mediate radii (r ~ 3 - 10 AU), behind the region shielded by the 
puffed-up inner rim. Gas temperatures at the inner rim are much 
higher when X-rays are included and, as a result, the inner rim 
is puffed up to higher and higher altitudes for increasing X-ray 
luminosities. As a direct extension of this theoretical work, the 
existence of the second bump could potentially be tested by con- 
tinuum imaging face-on protoplanetary disks in the near-infrared 
with, e.g., VLT or Keck. 
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Fig. 14. Radial column density distribution of H2O. FUV and X-ray fluxes are the same as Fig. [3] 



Scale height: Considering only FUV, we see that the scale 
height shows a maximum in the unattenuated parts (z/r > 0.5) 
of the disk. When X-rays are added we find that this maximum 
is smoothed over a larger region (out to r ~ 10 AU). The scale 
height in these regions is larger than one would expect from the 
flaring index in the outer regions of the disk. When moving to 
smaller relative height, e.g, z/r = 0.1, the break in the flaring 
index disappears. Another observational possibility would be to 
do continuum interferometry in the near-IR (VLTI) to directly 
measure the physical height of the inner rim. 



6.2. Chemical balance 

Ionization fraction: The ionization fraction reaches values as 
high as x e - ~ 10~ 2 in exceptional cases in our FUV-only mod- 
els, whereas X-rays easily maintain these ionization fractions 
throughout large portions of the disk. Even when the gas be- 
comes partially shielded, it can still maintain a significant ion- 
ization fraction, leading to an ion-molecule chemistry that can 
form molecules at low temperatures, which is not possible with 



neutral-neutral reactions as they usually have temperature barri- 
ers. 

Formation 0///2 through the H route: Formation of H2 dust 
is usually much more efficient in environments that have solar 
metallicities. Because of the high ionization fraction due to X- 
ray irradiation, the formation route H + H — > H2 + e is able 
to provide a significant addition of order 50% percent to the H2 
on dust formation route. Overall, the H2 to H abundance ratio is 
increased by at least two orders of magnitude when X-rays are 
present. 

Formation of water and OH: The OH and H2O abundances 
are more concentrated toward the inner regions of the disk, when 
only FUV is irradiating the disk. This is because the tempera- 
ture is only there sufficient to drive the neutral-neutral formation 
route. The outer disk shows significant enhancements (up to two 
orders of magnitude) when X-rays are added. The higher ion- 
ization fractions make it possible to form the species through 
ion-molecule reactions and sustain abundance levels of x H ,o ~ 
10~ 6 - 10~ 7 . Such abundance levels cannot be reached in outer 
disk models without X-rays. Only the warm inner disks allow 
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even higher levels of water abundance through warm neutral- 
neutral chemistry. 

Resulting abundance structures and radial column density 
profiles: Whereas neutral oxygen and CO are very stable to both 
FUV and X-rays, this is not the case for, e.g., C + , Ne + , and H2O. 
Ne + is strongly enhanced by X-rays and confined to the inner 
regions by larger FUV luminosities. This latter aspect will cer- 
tainly affect the line widths. A high L x /^fuv ratio is favorable 
for water formation, especially in the outer disk. 

6.3. Outlook 

In paper II, we will perform a radiation transfer analysis of the 
aforementioned species and correlate line fluxes and line widths 
directly to the FUV and X-ray luminosities. This will allow a 
discussion of the diagnostic value of these species and provide 
a theoretical framework for the interpretation of observational 
data. We will also discuss our results in the context of data 
obtained within several observational efforts (Spitzer, Herschel, 
and ground-based observing programs). 
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Appendix A: X-ray chemistry 

Our chemical network is composed of W s p=110 species (see Table |A. f| , linked 
by ~ 1500 reactions. 



Table A.l. Species included in the chemical network. 



1 10 species 



H, H + , H", H 2 , H+, H;, H+, He, He + , O, + , 2+ , 2 , OJ, 

OH, OH + , H 2 0, H z O + , H 3 + , CO, CO + ,C0 2 , CO+, HCO, HCO + , 

H 2 CO, 

N, N + , N 2+ , NO, NO + , C, C + , C 2+ , CH, CH + , CH 2 , CH+, CH 3 , CH+, 
CH 4 , CHJ, CH+, Si, Si + , Si 2+ , SiO, SiO + , SiH, SiH + , SiH+, SiOH + , 
S, S + , S 2+ , SO, SO + , S0 2 , SO+, OCS, CS, CS + , HS, HS + , HCS+, 
H 2 S + , H 3 S + , 

Mg, Mg + , Mg 2+ , Fe, Fe+, Fe 2+ , Ne, Ne + , Ne 2+ , Ar, Ar + , Ar 2+ , Na, 
Na + , Na 2+ , 

NH, NH+, NH 2 , NH+, N 2 H + , NH 3 , NH+, NH+, N 2 , HN+, CN, CN\ 
HCN, HCN+, HNC, C 2 H 2 , C 2 H+, HCNH + , CO#, B>0#, C0 2 #, 
CH 4 #, NH 3 #, 

PAH, PAH , PAH+, PAH 2+ , PAH 3 + 

References for the rate coefficients: Millar et al. ( 1986); Le nnon et al.| 
fT988} ; |Dffldini & Fossil l [TWT] ; |Anicich| ( |1993^|Badnell|j2006f ~ 



X-ray reaction rates, however, depend on the electron density n c , atomic hydro- 
gen density nu, and/or molecular hydrogen density nu 2 . These densities deter- 
mine how much of the absorbed X-ray photons will go into the different chan- 
nels for heating, ionization, and excitation. If the rate R shows a dependency 
on the particle densities, for example, the electron density (n e l), atomic hydro- 
gen density (nu), and/or molecular hydrogen density («h 2 ), the Jacobian can be 
expressed as: 



dFj 
drik 



dFj dFj dn c [ 
dn k <9«ei drib, 



dFj <5«h dFi dnu 2 
drm dn k <9«h? dnk 



(A.7) 



ProDiMo calculates this term analytically to ensure an accurate chemical solu- 
tion. Here we describe how the Jacobian terms for the X-ray reactions rates are 
calculated. We define the following quantities, which will be used below: 



tie] 



(A.8) 
(A.9) 



where nm\ is the total hydrogen nuclei density, Q(j,H) is the stoichiometric co- 
efficient for the hydrogen nuclei (which is for all those species that do not 
contain hydrogen, while for species containing hydrogen we have: Q(H,H)=l, 
Q(H 2 ,H)=2, Q(Hj ,H)=3, etc.). The electron density is denoted as n c \ and qj is 
the charge of the particle j. 



For a given species ;', the net formation rate reads (following Woitke et al. 2009): 



Jkl 



in jk 



Xsec' 
jk, 



(A.l) 



The terms involved are 



- Rjk-,n, the temperature-dependent rate for a two-body reaction where 
species i and / are formed out of species j and k 

- S?^ .j, a photo-reaction rate that depends on the local strength of the FUV 
radiation field 

- R^j k , a reaction that depends on the cosmic ray ionization rate 

- i? Xpr , and R 



the X-ray primary and secondary ionization reaction rates. 



Assuming statistical equilibrium = Oj, we obtain N sf non-linear equations 
for the N sp unknown particle densities n k . 



Fj(n k ) = 0. 



(A.2) 



For the species densities nk, this system of nonlinear equations is solved through 
a Newton-Raphson iterative method, which expands F, into a Taylor series in the 
neighborhood of 



ZdFj , 
— 6n k + 0(5n\) 



The code then needs to calculate the Jacobi matrix: 
dFt 



■'ik ■ 



drtk' 



(A.4) 



A.1. Primary ionization 

For the generic atomic species A, A + and molecular species AB, these are reac- 
tions of the kind: 



A + X ph - 


4 A 2+ 


A + +X ph - 


-> A 2+ 


AB + X ph - 


-> A 2+ 



+ 2e~ 
+ e~ 

+ B + 2e~ 
A + + B + + 2e" 
-> A + B 2+ + 2e" 



(A.10) 
(A.ll) 
(A. 12) 



where X p h is an X-ray photon. When the ('-th species is molecular ^A.12| , the 



Xnr 

reaction rate RV^ , ( for the reaction that destroys the species (' has two indices j, 
and /, for the resulting species. Otherwise the reactio n rate sho uld sim ply read 
R*^ T j as only one species is produced out of reactions , 



A.10 



and 



A.ll 



plicity, from now on, we use the molecular reaction rate as an example 
primary ionization. The rate is calculated as 



For sim- 
for X-ray 



„Xpr 

Vj1 



(E)F(E, r)dE [s" 



(A. 13) 



It depends on the X-ray radiation field F(E, r) at the point r where it is computed 
and on the cross section <x; of the species i. It does not depend on the local 
particle density. From equations |A.3| and |A.6| we see that the contribution to the 
equilibrium equation for the species i can be written as 



(A.3) Fj = F oth -Y J Rfl) r n i + Y J R > ; 



fl jl 

The contribution to the Jacobi element of this reaction is then 
dF olh 



Jii 



dn k 



(A. 14) 



(A.15) 



Neglecting terms of the order of 6n 2 k , the math SLATEC routines are used to find 
a set of fik that satisfies F,(nj + (5%) = 0: 

Jik ■ 5n k = -Fj. (A.5) 

When the rate of a given reaction does not depend on the particle density n k , the 
derivative of F, is straightforward, e.g, a photo-reaction that destroys species i: 



Jik ■ 



d_F\ 
dn k 



(A.6) 



A.2. Secondary ionization 

We consider X-ray secondary ionization for all the atomic species A and only 
for a single molecule H 2 : 



A + e 
H 2 +e~ 



A + + 2e~ 
H 2 + + 2e~ 
H + H + + 2e~ 



(A.16) 
(A.17) 
(A.18) 
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The reaction rate for the X-ray secondary ionization of species i is 



R X S ec = H X»(H) 
' Wn H 



[s-M. 



(A. 19) 



Substituting these results in Eq. |A.24| we obtain seven terms: 
dF Xxc „. 

-D K 



where r, is a parameter that takes into account the geometrical cross section 
of species i compared to the hydrogen cross secti on, ?1(h> is the total h ydrogen 
nuclei density, H\ is the X-ray energy deposition iMaloney et al. 1 996 1, and W 
[eV] is the mean energy consumed per ion pair I Dalga rnoetal.|1999) . 
W for hydrogen and all the other atomic elements except He is 



W H = W , H 1 + C, 



)a 



1 + C2 



«H 2 
I'H 



[eV], 



(A.20) 



where Wo_h is the collisional ionization rate for H in a pure neutral atomic gas 
(13.6 eV), M e i is the electron density, ;ih 2 is the molecular hydrogen density, and 
a, ci and ci are fitting parameters. 

Wh is then the energy needed to collisionally ionize hydrogen in a gas mix- 
ture with « e i , «h and «h 2 ■ The values for other at omic elements a re scaled up, 
considering the geometrical factor r, = tr™" / 'cr^ 11 Adamkovics et al. 2011 . W 
has been calculated for He and H2 as well; these cases win De treated separately 
further on. 



A3. Hydrogen 

The X-ray energy deposition is a quantity defined per unit of hydrogen nuclei as 
follows: 



H x = f 



<T tot (E)F(E,r)dE [erg^r 1 ], 



where <r tol is given by 



o-tot(£) = ^ 0";(£)«i7«<H> [cm 2 



(A.21) 



(A.22) 



and F(E, r) is the radiation field in unit of erg s" 1 cm" 2 eV~' . For a given atomic 
species (', at the equilibrium, the contribution of the secondary ionization to the 
volumetric rate is 



Hxnw H x n (H) 
= F oth -ri— m + rf— n,-. 



W H n H 



Following Eq. 



A.7 



the derivative of F ' c with respect to n k is 



tin/. 



(A.23) 



(A.24) 



m d(H x n (ll) ) tf x n (H >"/ <9W H H x n (W) nj H x n (i{) 

"H,* + „, Oik 



W H "h dn k W H n H Sn k Whi^ Wi«h 

Here we develop separately the partial derivative for all involved terms: 

■ Jz.jnjtTj(E)F(E,r)dE 

(r k (E)F(E, r)dE 
D K (A.25) 



d(H x n {H) ) _ Q 
dn k 



<9W H 
dn k 



d_ 
dii k 



W 0} 



H (i +C1 (^n(i + 



"<H>/ 

aq k n a - 1 I ,, H \ 

W).H ci -S— 1 + c 2 — 

<ID \ "H / 



«H 2 
" "H 



dnk 



1 j lJk 

W H n H 
H X n (H> 



(A.26) 



W H "h 



fym \"<H>/ I «h / 



Wj 

H x n (H) n 



W 2 n 2 \ V"<H)/ / m 



H x 

+ r,— — «/Wo.HCiQf, 1 
W^.n H \»<h>/ 



H 

Hx'HH)" 



"el V /. ^ «H 2 
I 1 1 + C 2 — I qt 



H x n {ii) n, 



— W ,H 1 + Ci 



"el 
"<H> 



C2(5h 2 



<5h,*- 



Wh" H 2 

The seven terms can be interpreted as follows: 

1 . The first term represents the variation of the energy deposition Hx when 
n k is increased. According to this term only, increasing n k brings more fast 
electrons in the gas phase via primary ionization of the species k. Hence 
more electrons are available for the secondary ionization of the species i. 

2. The second term directly impacts Rfj^f if i=k more particles are available 
for secondary ionization. 

3. The third calculates the variations in W if the species considered contains 
hydrogen. 

4. The fourth term comes from the dependency of Wh on «h (the higher «h 
the closer is Wh to Wo,h)- 

5. Analogous to the previous one, the fifth term shows that if the molecular hy- 
drogen density increases, atoms are less likely to be ionized by fast electrons 
because Wh will be higher. 

6. If the density of a positively charged particle increases, the electron density 
must increase at the same time as well. This favors Coulomb losses in en- 
ergy of the incoming fast electrons over secondary ionization of species i. If 
k =H~, this sixth term becomes negative because fewer electrons will be in 
the gas phase, favouring secondary ionization over Coulomb heating. 

7. The seventh term comes directly from the secondary ionization rate. The 
higher the atomic hydrogen density, the lower is Rfj^j ■ 

A.4. Molecular Hydrogen and Helium 

The mean energy per ion pair for H2 is 
Wh 2 = W ,h 2 (l+O (<,)") l + . 
while for helium, it is 

W He = W .He(l+Cl felf) 



(A.27) 



(A.28) 



where e* is 
el 



1.83&] 



1 + 0.83e d 



(A.29) 



with e c i is the electron fraction (n e l/"(H>)- Inserting these terms into |A.23| and 
generating the derivative with respect to n k gives the Jacobi elements tor molec- 
ular hydrogen and helium respectively. Since they look very similar to Eq. |A.26| 
we will not write them out. 



- Wo.hci 
+ W ,h 



(H) 

n" +l 
<H> 



1 + C2 - 



"H 



) 



"l c 2 
1 + c 'l 37- — S »2,k 



If 



W),H I 1 + C| I C 2 -^-<$Hjt- 
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Fig. B.l. Main heating sources throughout the disk: FUV luminosity increasing from Lfuy = 10 (left) to 10 erg/s. X-ray 
luminosity increasing from L x — to 10 32 erg/s. background heating by [CII] (blue), PAH heating (orange), photo-electric heating 
(dark purple), X-ray Coulomb heating (light red), heating by collisional de-exciation of H2 (blue-green), CI ionization heating 
(green-blue), infrared background by CO ro-vibrational lines (black), heating by thermal accomodation grains (white), cosmic ray 
heating (red), X-ray H2 dissociation heating (light blue), free-free absorption (green), background heating by Fell (light green), 
background heating by Sill (green-yellow), infrared background heating by H2O rotational transitions (yellow), heating by H2 
formation on dust (dark blue), and background heating by [OI] (purple). 

Appendix B: Additional figures 
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Fig. B.2. Main cooling sources throughout the disk: FUV luminosity increasing from Lft/v = 10 (left) to 10 erg/s. X-ray 
luminosity increasing from L x — to 10 32 erg/s. Lyman a cooling (black), [Fell] line cooling (blue-green), [OI] line cooling (white), 
[CII] line cooling (yellow), CO rotational and ro-vibrational cooling (red), H2O rotational cooling (green-blue), OH rotational 
cooling (light green), HCN line cooling (dark purple), [CI] line cooling (light purple), HNC line cooling (blue), cooling by thermal 
accomodation on grains (dark blue), CS line cooling (yellow-green), H2 line cooling (light blue), chemical cooling (red-orange), 
and free-free emission (green). 
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Fig. B.5. H2 abundances. FUV and X-ray fluxes are the same as Fig. 



B.3 
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Fig. B. 11. C + abundances. FUV and X-ray fluxes are the same as Fig. B.3 
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Fig. B.23. Ne abundances. FUV and X-ray fluxes are the same as Fig. B.3 
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